Localization-delocalization transition in disordered systems with a direction 
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Using the supersymmetry technique, we study the locahza- 
tion-delocahzation transition in quasi-one-dimensional non- 
Hermitian systems with a direction. In contrast to chains, 
our model captures the diffusive character of carriers' mo- 
tion at short distances. We calculate the joint probability 
of complex eigenvalues and some other correlation functions. 
We find that the transition is abrupt and occurs as a result 
of an interplay between two saddle-points in the free energy 
functional. 

72.15.Rn, 73.20.Fz, 74.60.Ge 



Non-Hermitian models with disorder have attracted re- 
cently a considerable attention p|-pl]|. Non-Hermitian 
Hamiltonians appear in context of flux lines in supercon- 
ductors j^], in transfer phenomena in lossy media in 
hydrodynamics 0], in QCD |l0| , in quantum mechanics 
of open systems . 

The most important property of the non-Hermitian 
Hamiltonians is that their eigenenergies can be complex 
independently of their origin. Very interesting physical 
systems are models with a direction. The simplest Hamil- 
tonian TL of such models is written as 



„^(P+ih):^ 

2 m 



(1) 



where p is the momentum operator, V (r) is a random 
potential and h is a constant vector. One comes to the 
Hamiltonian 7i, Eq. (|^), reducing the d+1 dimensional 
classical problem of vortices in a superconductor with 
columnar defects to a d-dimensional quantum problem. 
Then, the vector h is proportional to the component of 
the magnetic field perpendicular to the line defects Q]. 

Considering the one-dimensional (ID) version of the 
Hamiltonian H, Eq. (P, the authors of Ref. [Q pre- 
dicted a localization-delocalization transition that occurs 
at he = , where Ic is the localization length at h =0. 
At h < the "imaginary vector-potential" h can be 
removed by the "gauge transformation" 

</>/c (r) = (/>ofc (r) exp (hr) , 4>k (r) = (j)Qk (r) exp (-hr) 

(2) 

In Eq. (|^), (pk (r) and 0^ (r) are the right and left eigen- 
functions, 0ofe (r) are those at h =0. Carrying out the 
transformation, Eq. (^, we see that the eigenvalues of 
the localized states eu do not depend on h. However, the 



gauge transformation applies only if the functions (pk (r) 
do not grow at infinity. For eigenfunctions of the form 
00 ~ exp(— |r|/^c), this is possible ii h < he- For h > he, 
one should take extended wave functions, which leads to 
complex eigenenergies. 

Qualitative arguments of Ref. |l[ have been confirmed 
for the disordered chains numerically as well as analyti- 
cally [0-^ . At the same time, they are quite general and 
one might expect that such a transition occurs in more 
complex systems provided all states at h =0 are local- 
ized. For example, one can think of disordered quasi- 
one-dimensional wires and films. Using the mapping of 
Ref. one expects such models to be relevant for de- 
scription of the vortices in slabs and 3D bulk supercon- 
ductors. The wires and films are richer than the chains 
because at distances larger than the mean free path I but 
smaller than the localization length Lc the motion is dif- 
fusive, whereas in disordered chains only ballistic motion 
or localization are possible. As a result, properties of the 
localization-delocalization transition in wires and films 
can differ. The problem has not been addressed yet and 
our paper is aimed to provide its basic understanding. 

We study the distribution function of complex eigenen- 
ergies P (e, y), where e and y are the real and imaginary 
parts of the energy, and correlations of eigenfunctions of 
the Hamiltonian Ti., Eq. (|^), at different points for wires 
using the supersymmetry technique p^ . A proper non- 
linear supermatrix cr-model has been derived in Ref. 
Its free energy functional Fh [Q] reads 



Fh [Q] = Y^^^ I - 4 (7A + yAiTa) Q]dr 



(3) 



where D is the classical diffusion coefficient, 1/ is the den- 
sity of states, and the standard notations for the super- 
trace Str and supermatrices Q, T3, A, and Ai are used 
jl|,|. The free energy functional Fh [Q], Eq. @, differs 
from the conventional one Fq [Q] (describing the Hermi- 
tian problem) by the presence of the "gauge-invariant" 
combination dQ = VQ + h [Q, A] instead of VQ. 

The zero-dimensional (OD) version of the cr-model is 
equivalent to models of random non-Hermitian matrices 
Q , where no transition occurs. But the cr-model, Eq. (|^), 
holds in any dimension and one might expect that already 
the ID version would manifest the transition. So, one 
could describe the transition using the transfer-matrix 
technique developed for the cr-models p^ , p^ . 

Surprisingly, the ID cr-model with Fh [Q], Eq. (|^), does 



1 



not have any transition and physical quantities of interest 
are smooth functions of h giving always a finite proba- 
bility of complex eigenenergies. At first glance, this re- 
sult is in an evident contradiction with the arguments 
of Ref. The resolution of this paradox is quite in- 
teresting. It turns out that the ct- model with 
Eq. (||) , describes properly only the region of the delocal- 
ized states. The functional Fo[Q] should be used in the 
localized regime, for any finite h < he, where he is a criti- 
cal vector-potential. This replacement of the free energy, 
resembling a first order transitions, results in an abrupt 
change of the distribution of eigenvalues. At h < he, all 
eigenenergies are real whereas at h > he one gets a broad 
distribution in the complex plane. In the limit h ^ he, 
our results are universal for any dimension. All this is in 
a contrast with the results obtained for the chains, where 
the spectrum changes smoothly, showing the appearance 
of "arcs" when increasing h |]l],|j , or "wings" when start- 
ing from the regime of the strong non-Hermiticity j^]. 

Having formulated the rough picture of what happens 
when changing the parameter h, let us present some de- 
tails of calculations. We first introduce the quantities 
which are to be studied. Since eigenvalues of the Hamil- 
tonian (|]) can be complex, it is convenient to double the 
size of relevant matrices |^ , thus "hermitizing" the prob- 
lem. In such an approach, the role of the Green function 
is played by the function 



70fe(r)0jr') 



(4) 



with an eigenenergy ek = ej^ + i^'^- The "density-density 
correlator" in the present context is given by 

y(r, r'; e, y) - - lim (i?(r, r')B{v', r)) , (5) 

IT 7^0 

where brackets imply impurity averaging. The limit 
7 — > 0, understood in all correlators, becomes important 
as soon as /i — > 0. The function Y{r,r';e,y), Eq. (||), 
establishes a link between the localization properties and 
the joint probability density of complex eigenenergies 



Pie,y)^^ J drdr'r(r,r';e,y) 



(6) 



where V is the volume. We introduce also another impor- 
tant correlator X (r, r'; e,y) — C (r, r'; e, y)+C (r', r;e,y) , 

C (r, r'; e,y) ^ ^ lim {B{v, r')B*{r, r')) (7) 

ZTT 7^0 

The correlator Y{r,r'; e,y), Eq. is invariant under 
the transformation, Eq. (J^, but X (r, r'; e, y) is not. 

We further express as usual |^,^ the correlation func- 
tions in terms of integrals over eight-component super- 
field ip{r), average over the white- noise disorder poten- 
tial V{r), decouple -0^ term by 8 x 8 matrix Q{r) and 



integrate over ip (r). As a result, we obtain a functional 
integral over Q with a free energy functional 



F[Q] = \j ST 



r I ; In — «7t 

4r 



Q(r) 
2t 



dr 



n 



(8) 



where Ti.' = — ihV/m + yr^, and r is the mean free time. 

The next standard step is to find the minimum of F[Q] 
neglecting Ti.' . The minimum is reached at 



= VAV 



(9) 



VV ~ 1, y+ = KVK. (The notations are the same as 
in Ref. |]3|,|l2|). Expanding the functional F[Q] near the 
minimum in the gradients of Q and in Ti.' one comes to 
the functional Fh [Q], Eq. . This is exactly the way how 
the cr- model, Eq. was derived in Ref. However, 
it has been mentioned that this cr-model is not valid in 
the localized regime. What is wrong in the derivation? 

It not difficult to find that the functional F [Q] has also 
another minimum. Performing the transformation 



Q = exp (Airh) Q exp (— Airh) 



(10) 



we rewrite the functional F [Q] , Eq. (||) in terms of Q. 
As a result, the imaginary vector-potential h is removed 
from F and the minimum is achieved at 



Q = VAV 



(11) 



which corresponds to Q varying in space. 

Which of these two saddle points should be cho- 
sen? The answer depends on the value of /i. To clar- 
ify this question we consider the correlation functions 
Y (r, r'; e, y) and X (r, r'; e, y), Eqs. (||J^). They can be 
written as functional integrals over the supermatrices Q 



(12) 



where the sign + (— ) corresponds to the correlator Y 
{X), gi± = Q"(r)±Q''(r), = Qi2(i.)±g2i(i.), and 
the matrices Q' are taken at r'. The symbol {■■)q stands 
for averaging with the functional F [Q] , Eq. (||) . To cal- 
culate the integral in Eq. ( p2| ) let us use the transforma- 
tion, Eq. (pO[), and take the saddle-point, Eq. (pT]). Since 
the combination of the supermatrices Q entering Eq. ( p^ ) 
for the function Y is invariant under the transformation, 
Eq. ([l^), h is gauged out in this function. Therefore, one 
can use the standard results of the transfer-matrix ap- 
proach developed for the Hermitian case. The final 
result for the correlator Y (r, r'; e, y) 



Y {r,v';e,y) ^ Poo (r) S (y) 



(13) 



where r = |r — r'|, contains the function poc (r) charac- 
terizing localization properties 
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Poo (r) 



El 

A: 



bok{r)f\<l^ok{r')fS{e^e,) (14) 



In disordered wires, the function poo (t) describes the de- 
cay of the wave functions at r ^ Lc 



1 



Po 



4L, 
r 



3/2 



exp 



r 



(15) 



where Lc is the locahzation length. For the unitary en- 
semble it equals Lc = 2t:i'SD, S is the cross-section. 

In contrast, the vector-potential h enters explicitly the 
pre-exponential of the function X after making the trans- 
formation, Eq. (p^. However, the dependence on h is 
simple and calculations for X (r, r'; e, y) can be performed 
in the same way as for Y (r, r'; e, y) yielding 



X (r, r'; e, y) = cosh {2hr) poo (r) S (y) 



(16) 



Eqs. ( |l^ , [T^ ) demonstrate that, provided one may per- 
form the transformation, Eq. ([To|), all eigenenergies re- 
main real. However, the validity of this procedure de- 
pends on the value of h. The function X (r, r'; e, y) does 
not grow at infinity only \i h < where 



(8L,)-i 



(17) 



At h > he Eqs. (^,^6|) cannot be used because this would 
correspond to growing wave functions, which are forbid- 
den for a closed geometry. This agrees with the argu- 
ments of Ref . ||l| . According to Ref . one should use at 
h > he extended states having complex eigenenergies. 

In the present formalism, the other saddle-point, 
Eq. (H), of the free energy F[Q] should be taken in the 
regime h > he- This leads to the cr-model in the form of 
Eq. (H). Expecting that the eigenenergies become com- 
plex, we determine their distribution function P (e, y), 
Eq. (H). Following the transfer-matrix technique [ p^JT2| 
we write the function P (e, y) in the form 



Pie,y) 



* (Q) {QltP24 ~ QltP2A) dQ (18) 



In Eq. ([l^), ^ (Q) is the partition function of the semi- 
infinite wire, the matrix function P is the partition func- 
tion between the points r and r' multiplied by Q''i'{Q') 
and integrated over r' . As usual, proper differential equa- 
tions for and P are derived from F[Q], Eq. @. 

In order to carry out these calculations it is necessary 
to choose a parametrization of the supermatrices Q. For 
the non-Hermitian problem involved the parametrization 
of Ref. 1^ is most natural. For simplicity we consider 
the unitary ensemble, where the supermatrices Q can be 
parametrized in the form 



Q^T(0]QoT 



Qq 



cos (p 
-T^ sin (fi 



-T3 sm (fi 
— cos (fi 



f = 







§ = 


\ 








[oieij 



(19) 



In Eqs. (|19[), the supermatrices T{0) contain not only real 
variables 9 but also Grassmann variables. The "angles" 
vary in the following intervals — 7r/2 < ip < tt/2, ~oo < 
X < 00, — TT < d < TT, —00 < 9i < 00. 

The variables and (p are not equivalent. For example, 
neglecting the gradient terms in Eq. (||) one comes to the 
OD free energy F^^\ containing the variables only 

P(°)[(y3] = h^{Xi - ty/2hY + + yl2h'^ f , (20) 

where Ai = sinhx, A — sin(y5, h^ = 2ft,^L^ and y = 
2yL^/D. In contrast to real magnetic field H, gradually 
suppressing and finally freezing out some degrees of free- 
dom with increasing H, the imaginary magnetic field h 
shifts the saddle point as a whole. Noticeable changes in 
behavior occur only at zLy > 2h^, where siiap « 1. 

Calculations performed in Ref. |^ for the OD case show 
that the variables 9 play a minor part in OD. They do not 
enter P^"-* [ip] but their role is even less pronounced due 
to the singularity of the Jacobian at 9, 0i — > 0. A de- 
tailed discussion of Ref. § leads to the conclusion that 
one should replace the ^-dependent part of the Jacobian 
by a constant and put everywhere else 9 — 0. It is also 
interesting to note that the free energy Fh \Q], Eq. (||), 



is not invariant against the replacement, Eq. (10). Using 
the parametrization, Eqs. (^9|), one obtains that this re- 
placement leads to the shift 9 = 9 — 2irh, 9i = 9i — 2rh. 
This shift changes the contour of the integration over 
— TT < < TT in a complicated way, thus demonstrating 
the violation of the "gauge symmetry" . 

Effective disappearance of the variables 9 shows that 
the correspondence between the problem under consider- 
ation and random flux models [|6| is rather obscure be- 
cause in the latter, the variables 9 play a crucial role. For 
the quantities we have studied, Eqs i,|,|), this mapping 
is inadequate. 

As in the OD case, one should put everywhere 6, 61 = 
when deriving the transfer matrix equations. As a result, 
differential equations for '5 and P contain A and Ai only 

7^^» = i^^W^r, 'HP+ = Lc{iXi - iX)-^ , (21) 

n - ^dxil - X')Jodx + ^dx,{l + Xl)Jodx, , 
•Ja Jo 

where Jq = l/(Ai + iX)'^. 

Analysis of Eqs. ( pT| ) for hLc ^ 1 is to some extent 
similar to the one in the limit of high frequencies for the 
conventional problem of localization |13] . We obtain 



^' ~ exp[-h{Xi+iy/2h^f -h{X-y/2hy 
which gives after substitution into Eq. (11^ 



P^ 



1, 
0, 



1^1 < 2h^ 

\y\ > 



'5/2/1 



(22) 
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The form of the functfon P(e, y), Eq. (p^, is the same as 
the OD result of Refs. This result does not depend 

on the dimensionality and corresponds to the elliptic law 
for strongly non-Hermitian random matrices JlSf . 
Analogous calculations for Y, Eqs. (Q, ^) yield 

Y{r) ~ i//3exp(-2/3|r|) , /? = y^{h^ - y^/ADK^) . (23) 

Analytical study of Eqs. is hardly feasible at 

hLc ~ 1. To solve them numerically, we use the stan- 
dard over-relaxation method with Chebyshev accelera- 
tion. The results for the distribution P(e, y) are pre- 
sented in Fig. 1 for several values of h. The lowest value 
oi h — l/4-\/2 corresponds to the critical /ic, Eq. (|l7[). 
For comparison we present the OD result of Refs. [||J^ 

The curves for the ID and OD cases are rather close to 
each other. However, their evolution with decreasing h 
is drastically different. The OD curve tends smoothly to 
the (^-function v5 (y) when h 0, wheareas the ID curve 
changes abruptly to this expression at h = he- In the re- 
gion h > he the states are extended and the exponential 
decay of the correlation function Y (r), Eq. (p3|), is ob- 
tained after summation over many states with different 
phase differences. 

FIG. 1. The joint distribution function of complex eigen- 
values P{e), Eqs. for the quasi- ID infinite wire (solid line) 
and the OD case (dashed fine), for h = /ic, 0.2, 0.25. 

P(e, y) I ■ ■ ■ ■ ■ ■ 1 




In conclusion, we studied the localization-delocaliza- 
tion transition for non-Hermitian disordered systems 
with a direction, Eq. (Q) and broken time-reversal sym- 
metry. We found that the transition occurs due to an 
interplay between two saddle points in the free energy 
_F[(5], which results in the abrupt transition at the criti- 
cal field h = he. Below this field, all eigenvalues are real 
and wave functions are localized. Above the critical field 
the eigenvalues form a broad distribution in the complex 
plane. We believe that the abrupt transition found is 
quite natural for the original problem of the vortex lines 
in the presence of columnar defects j^]. Provided the ge- 
ometry of the sample is closed and it is infinitely long, 



vortex lines should abruptly change their behavior, from 
being pinned by columns to become oriented along the 
magnetic field. 
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